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Abstract. 

The Tartarus database contains products for 529 ASCA observations of active galactic nuclei. We have been updating 
Tartarus to include observing sequences conducted after 1999 January. The revised database will contain products for 375 
objects, with a total of 614 observing sequences. We have begun a systematic timing analysis of the Tartarus light curves. 
We present here some preliminary results of an investigation into the relation between excess variance and black-hole mass. 
Having optimised our analysis to minimize the scatter in the variance measurements, we find that the narrow-line active 
galactic nuclei follow roughly the same relation as the broad-line objects. 



X-RAY VARIABILITY IN ACTIVE 
GALACTIC NUCLEI 

The existence of variability in the X-ray emission from 
active galactic nuclei (AGNs), on time-scales of days to 
years, was established roughly three-decades ago [e.g. 1, 
and references therein]. Variations on time-scales shorter 
than this are also common. The launch of EXOSAT, with 
its long orbit, meant that well-sampled X-ray light curves 
could be obtained for time-scales of minues to days. 
The power-spectra generated from EXOSAT light curves 
could be described by an unbroken power-law [e.g. 2, 3]. 
Archival observations spanning weeks to months showed 
that, for two objects, the power-law flattened at low fre- 
quencies [4, 5]. More recently, RXTE and XMM-Newton 
observations have shown that the power-spectra of AGN 
light curves exhibit broken power-law shapes very sim- 
ilar to those seen in galactic black-hole X-ray binaries 
[e.g. 6, 7]. 

As well as power-spectral density, the excess variance 
can be used to characterise the X-ray variability [e.g. 8]. 
Excess variance is defined as the variance in the light 
curve not due to statistical fluctuations. The window 
function affects the excess variance less than it does the 
power-spectral density. Therefore, the excess variance is 
an especially useful quantity for the description of light 
curves containing many gaps. 



An anti-correlation was found between excess vari- 
ance and luminosity for a sample of AGNs observed by 
the Advanced Satellite for Cosmology and Astrophysics 
(ASCA) [8]. Later work using ASCA data found that, for 
a given luminosity, the X-ray light curves of narrow-line 
Seyfert 1 (NLS1) galaxies exhibit a larger excess vari- 
ance than broad-line Seyfert 1 (BLS1) galaxies [9, 10]. 
This result can be explained if, for a given luminosity, 
the NLS 1 s have a lower black-hole mass than the BLS Is. 
A lower mass results in faster variations, and this is seen 
as a larger excess variance. Alternatively, other proper- 
ties of the variability, not just the time-scale, might differ 
between NLS Is and BLS Is. These alternatives can be 
tested through examining the dependence of the excess 
variance on black-hole mass. 

The relationship between excess variance and black- 
hole mass, for Seyfert 1 galaxies and quasi-stellar objects 
(QSOs), has been studied using ASCA data (time-scale ~ 
1 d) [11, 12]. These studies revealed an anti-correlation 
between mass and variance. Moreover, the narrow-line 
(NL AGNs) and broad-line (BL AGNs) objects followed 
the same relation. This relation could be explained if the 
power-spectrum is an unbroken power-law [P(v)°< V~ a ] 
with a ~ 2 [11]. 

An investigation using RXTE data (time-scale ~ 
300 d) presents a different picture [13]. The BLS Is, and 
the NLS1 MCG-6-30-15, follow a variance-mass rela- 



tion that is consistent with a universal power-spectral 
shape. This 'average' spectrum comprises a twice- 
broken power-law which breaks from a = to 1 at the 
so-called iow-frequency break', and from a = 1 to 2 
at the 'high-frequency break' . The break frequencies are 
inversely proportional to black-hole mass, and the am- 
plitude, when represented in power x frequency space, is 
constant. This power-spectral shape resembles that seen 
in the black-hole binary Cyg X-l in the low/hard state. 
The other NLS1 in the sample, NGC 4051, did not fol- 
low the above relation. The excess variance of NGC 405 1 
was consistent with a singly-broken power-law, breaking 
from a = 1 to 2, with a break frequency that is 20 times 
higher than for BLSls. This shape is similar to that seen 
in Cyg X-l in the high/soft state. A power-spectral anal- 
ysis of NGC 4051, using RXTE (time-scale ~ 6.5 years) 
and XMM-Newton (time-scale ^100 ks) data, confirms 
that NGC 4051 is indeed an analogue to high state bi- 
naries [6]. Moreover, in plots of high-frequency break 
versus black-hole mass, the NLSls were found to have 
systematically higher break frequencies than the BLSls 
[13,6]. 

The excess variance versus luminosity relation for 
the BLS1 galaxies can be explained as a consequence 
of the variance-mass relation if the BLSls radiate at 
~ 13 per cent of the Eddington luminosity [13]. The 
location of NGC 405 1 in the variance-luminosity plane 
is consistent with that object radiating at ~ 10 per cent 
Eddington [13]. 

These results from RXTE and XMM-Newton suggest 
that the power-spectral properties of AGN depend on 
black-hole mass and at least one other parameter. The de- 
pendence with mass is perhaps the least surprising, since 
we expect the size of the emission region — and, there- 
fore, the variability time-scale — to scale with mass. A 
useful step forward, then, is to account for this mass de- 
pendence and investigate how the X-ray variability de- 
pends on other parameters. This is best achieved through 
maximising the sample size and minimising observa- 
tional scatter. 



TARTARUS: A DATABASE OF ASCA 
OBSERVATIONS OF AGN 

The ASCA satellite was launched into a 96 minute or- 
bit on 1993 February 20 and ceased observing on 2000 
July 14 [14] 1 . The satellite comprised four X-ray tele- 
scopes and four corresponding focal-plane instruments 
[15]. There were two solid-state imaging spectrometers 



(SIS0 and SIS 1) [16] and two gas imaging spectrome- 
ters (GIS2 and GIS3) [17]. The effective energy ranges 
were 0.5-10 keV and 0.7-10 keV for the SIS and GIS 
instruments, respectively. 

The Tartarus database contains products for all ASCA 
observations with targets designated as AGN [18]. The 
products comprise SIS and GIS event files, images, spec- 
tra, and light curves. Version 2 is presently available 2 at 
Goddard Space Flight Centre (GSFC), and consists of 
529 observing sequences. We are updating and improv- 
ing the database to include data collected after 1999 Jan- 
uary. The revised database will consist of 614 sequences 
with a total of 375 objects. Version 3 of Tartarus is to be 
made available at GSFC and Imperial College London. 

We are conducting a systematic timing analysis of the 
light curves in the Tartarus database. Tartarus is well 
suited to such a project because it contains data for a 
large sample of objects. The overall aim of the project 
is to uncover any correlations that exist between the 
properties of the X-ray variability and other parameters. 

We have begun the project by investigating those AGN 
for which the black-hole mass has been estimated. We 
present here some preliminary results of an investigation 
between excess variance and black-hole mass. Our work 
is similar to that conducted using RXTE data as described 
in the previous section. 



ANALYSIS 

We restricted our analysis to narrow- and broad-line 
AGN. Our sample contains Seyfert 1 galaxies and QSOs. 
We do not differentiate between Seyfert galaxies and 
QSOs, but we shall simply refer to each object as ei- 
ther a narrow- or broad-line AGN [12]. We used the 'rms 
spectrum' [24] reverberation-mapping mass estimate if 
available, otherwise we used the mass estimate deter- 
mined from the optical luminosity [e.g 23]. The objects 
in our sample, along with their black-hole masses, are 
listed in Table 1 . The mass measurements with no error 
bars were estimated via the optical luminosity. The un- 
certainty in these measurements is about 0.5 dex [12]. 
Objects with Hj3 < 2000 km s~ l were classified as NL 
AGNs [9, 12, 25,26]. 

The variance in a light curve is equal to the integral of 
the power-spectral density between the frequency limits 
set by the time resolution and duration of the light curve. 
Therefore, the use of different binnings and durations per 
object will introduce artificial scatter. We wish to mini- 
mize this scatter. Therefore, we used a time resolution of 
256 s for all light curves, and restricted them to be be- 



see also http://heasarc.gsfc.nasa.gov/docs/asca/ 
ascagof . html 



2 http// : tartarus . gsf c . nasa . gov 



TABLE 1. Name, type and black-hole mass for each ob- 
ject having at least 1 valid Tartarus light curve. The objects 
exhibiting significant variability are listed first. 



Object 


Type 


Mass 


Ref. 


Name 
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NGC4051 


NL AGN 


05+ 06 


[19] 


MRK 766 


NL AGN 


0.0832 


[12] 


AKN 564 


NL AGN 


0.115 


[12] 


MCG-6-30-15 


NL AGN 


0.155 


[12] 


MRK 1040 


BL AGN 
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NGC 3783 


BL AGN 


0.87 + 0.11 


T221 


MRK 290 


BL AGN 


1.12 


[12] 


NGC 4151 


BL AGN 


1.20+^3 

—0.70 


[21] 


TON SI 80 


NL AGN 


1.23 


[12] 


NGC 3516 


BL AGN 


1.68 + 0.33 


[20] 


NGC 3227 


BL AGN 


3.6+1.4 


[20] 


MRK 509 


BL AGN 


9.2+1.1 


[21] 


NGC 5548 


BL AGN 
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PG 0026+12 


NL AGN 


2 - 66 -0.55 


[21] 


MRK 279 


BL AGN 


4.17 


[12] 


F9 


BL AGN 


8.3+^ 


[21] 


MRK 841 


BL AGN 


14.1 


[12] 


PG 1116+215 


BL AGN 


16.2 


[23] 


PG 0804+761 


BL AGN 


16.3+1-g 


[21] 


AKN 120 


BL AGN 
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[21] 


MCG-2-58-22 


BL AGN 


34.7 


[12] 


MRK 1383 


BL AGN 


37 +13 
-"-16 


[21] 


IRAS F04250-5718 


BL AGN 


38.0 


[12] 


HE 1029-1401 


BL AGN 


120 


[23] 



tween 30 and 40 ks in duration. 

A sample of measurements of the excess variance of a 
red-noise process will possess a large scatter about their 
mean [27, see also the following section]. This scatter 
is, of course, reduced by averaging over many realiza- 
tions. Therefore, in our analysis, light curves longer than 
40 ks were subdivided into many shorter segments. This 
allowed us to maximize the number of variance measure- 
ments per object. 

We extracted a 2-10 keV combined SIS0 and SIS1 
light curve for each object. We only used light curves 
comprising at least 20 bins, with each bin having no 
fewer than 20 counts. Some light curves, having being 
initially extracted, contained some bins having fewer 
than 20 counts. In these cases, we excised bins from the 
light curve, on the basis of fractional exposure, so that 
only valid bins remained. The invalid bins could thus be 
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log (Block Hole Moss / M e ) 

FIGURE 1. Logio of excess variance versus logio of black- 
hole mass. Open and filled circles represent NL AGNs and 
BL AGNs, respectively. The stars indicate the masses of those 
objects in which variability was not detected. The solid line 
shows the expected variance assuming a twice-broken power 
law, and the dashed lines show the range in the expected value. 



removed from the light curve without introducing bias. 

There were 29 objects for which we had at least a sin- 
gle valid light curve segment. These objects are listed in 
Table 1 . We determined the % 2 corresponding to the hy- 
pothesis of a constant counting rate for each segment. We 
then summed all of the j 2 s and degrees-of-freedom for 
each object, and thereby tested whether that object ex- 
hibited variability. We detected variability in 16 objects 
at the 95 per cent confidence significance level (see Ta- 
ble 1). We are yet to determine the upper limits for those 
objects in which variability was not detected. 

For the 8 objects having at least 5 measurements of ex- 
cess variance, the uncertainty in the mean variance was 
measured from the observed scatter in those measure- 
ments. The standard deviation of the scatter about the 
mean was found to be <~ 60 per cent of the mean. We 
used this relation to assign uncertainties for the 8 objects 
having fewer than 5 variance measurements. 



RESULTS, DISCUSSION AND FUTURE 
DIRECTIONS 

A plot of logio of the excess variance versus logio 
of the black-hole mass is presented in Fig. 1 . The open 
and filled circles indicate the NL AGNs and BL AGNs, 
respectively. The stars indicate the masses of the 13 
objects in we did not detect variability. 

The solid line in Fig. 1 shows the expected variance 
assuming a twice-broken power-law. We determined this 
using the relation found from RXTE data for BLSls [13], 



except that we also accounted for red-noise leak. The 
dashed lines indicate the rough range of the expected 
variance based on the errors from the RXTE analysis. 

The NL AGNs and BL AGNs follow roughly the same 
relation. This conclusion was also reached with previous 
analyses of ASCA data [11, 12]. We note again that 
our analysis is an improvement over the previous work 
because we have calculated the excess variance from 
light curves having nearly the same length. 

The next step in our work is to include the GIS data in 
our analysis to improve our signal-to-noise and increase 
the number of valid light curves. We will also investigate 
the variance-mass relationship in different energy bands. 

Most importantly, our future work will involve simu- 
lations similar to those used in modeling AGN power- 
spectra [e.g. 7]. This will allow us to account for the 
sampling of the light curves and, thereby, permit a robust 
comparison with a hypothesised universal power spec- 
trum. Moreover, we shall quantify the extent to which 
the mean excess variance of each object differs from the 
expected variance. These residuals can then be examined 
against other parameters, such as mass-accretion rate. In 
this way, we can investigate the extent to which various 
parameters, other than black-hole mass, are associated 
with the X-ray variability. 



ACKNOWLEDGMENTS 

PMO acknowledges financial support from PPARC. This 
research has made use of the TARTARUS database, 
which is supported by Jane Turner and Kirpal Nandra 
under NASA grants NAG5-7385 and NAG5-7067. 



REFERENCES 



10. Leighly, K. M, ApJS, 125, 297-316 (1999). 

11. Lu, Y, and Yu, Q., MNRAS, 324, 653-658 (2001). 

12. Bian, W., and Zhao, Y, MNRAS, 343, 164-168 (2003). 

13. Papadakis, I. E„ MNRAS, 348, 207-213 (2004). 

14. Tanaka, Y, Inoue, H., and Holt, S. S., PASJ, 46, L37-L41 
(1994). 

15. Serlemitsos, P., Jalota, L., Soong, Y, Kunieda, H., Tawara, 
Y, Tsusaka, Y, Suzuki, H., Sakima, Y, Yamazaki, T., 
Yoshioka, H., Furuzawa, A., Yamashita, K., Awaki, H., 
Itoh, M., Ogasaka, Y, Honda, H., and Uchibori, Y, PASJ, 
47, 105-114 (1995). 

16. Gendreau, K., X-ray CCDs for Space Applications: 
Calibration, Radiation Hardness, and Use for Measuring 
the Spectrum of the Cosmic X-ray Background, Ph.D. 
thesis, Massachusetts Institute of Technology (1995). 

17. Ohashi, T., Ebisawa, K., Fukazawa, Y, Hiyoshi, K., Horii, 
M., Ikebe, Y, Ikeda, H., Inoue, H, Ishida, M., Ishisaki, 
Y, Ishizuka, T., Kamijo, S., Kaneda, H., Kohmura, Y, 
Makishima, K., Mihara, T., Tashiro, M., Murakami, T., 
Shoumura, R., Tanaka, Y, Ueda, Y, Taguchi, K., Tsuru, 
T, and Takeshima, T, PASJ, 48, 157-170 (1996). 

18. Turner, T. J., Nandra, K., Turcan, D., and George, I. M., 
"Tartarus- An ASCA AGN database," in AIP Conf. Proc. 
599: X-ray Astronomy: Stellar Endpoints, AGN, and the 
Diffuse X-ray Background, 2001, pp. 991-994. 

19. Shemmer, O., Uttley, P., Netzer, H., and McHardy, I. M., 
MNRAS, 343, 1341-1347 (2003). 

20. Onken, C. A., Peterson, B. M., Dietrich, M., Robinson, 

A. , and Salamanca, I. M., ApJ, 585, 121-127 (2003). 

21. Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, 

B. T, and Giveon, U., ApJ, 533, 631-649 (2000). 

22. Onken, C. A., and Peterson, B. M., ApJ, 572, 746-752 
(2003). 

23. Woo, J. H., and Urry, C. M., ApJ, 579, 530-544 (2002). 

24. Peterson, B. M., Wanders, I., Bertram, R., Hunley, J. F., 
Pogge, R. W., and Wagner, R. M., ApJ, 501, 82-93 (1998). 

25. Boroson, T. A., and Green, R. F., ApJS, 80, 109-135 
(1992). 

26. McLure, R. J., and Dunlop, J. S., MNRAS, 327, 199-207 
(2001). 

27. Vaughan, S., Edelson, R., Warwick, R. S., and Uttley, P., 
MNRAS, 345, 1271-1284 (2003). 



1 . Marshall, N., Warwick, R. S., and Pounds, K. A., MNRAS, 
194, 987-1002 (1981). 

2. Lawrence, A., and Papadakis, I., ApJ, 414, L85-L88 
(1993). 

3. Green, A. R., McHardy, I. M., and Lehto, H. J., MNRAS, 
265, 664-680 (1993). 

4. McHardy, I., Mem. It. Astr. Soc, 59, 239-259 (1988). 

5. Papadakis, I. E„ and McHardy, I. M., MNRAS, 273, 
923-939 (1995). 

6. McHardy, I. M., Papadakis, I. E., Uttley, P., Page, 
M. J., and Mason, K. O., MNRAS, accepted (2003), 
[astro-ph/0311220]. 

7. Markowitz, A., Edelson, R., Vaughan, S., Uttley, P., 
George, I. M., Griffiths, R. E., Kaspi, S., Lawrence, A., 
McHardy, I., Nandra, K., Pounds, K., Reeves, J., Schurch, 
N, and Warwick, R., ApJ, 593, 96-114 (2003). 

8. Nandra, K., George, I. M., Mushotzky, R. F, Turner, T. J., 
and Yaqoob, T, ApJ, 476, 70-82 (1997). 

9. Turner, T. J., George, I. M., Nandra, K., and Turcan, D., 
ApJ, 524, 667-673 (1999). 



